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ABSTRACT 

We apply the intrinsically symmetrical, decelerating relativistic jet model developed by 
Laing & Bridle for 3C 31 to deep, full-synthesis 8.4-GHz VLA imaging of the two low- 
luminosity radio galaxies B2 0326+39 and B2 1553+24. After some modifications to 
the functional forms used to describe the geometry, velocity, emissivity and magnetic- 
field structure, these models can accurately fit our data in both total intensity and 
linear polarization. We conclude that the jets in B2 0326+39 and B2 1553+24 are at 
angles of 64° ± 5° and 7.7° ± 1.3° to the line of sight, respectively. In both objects, we 
find that the jets decelerate from 0.7 - 0.8c to <0.2c over a distance of approximately 
lOkpc, although in B2 1553+24 this transition occurs much further from the nucleus 
than in B2 0326+39 or 3C31. The longitudinal emissivity profiles can be divided into 
sections, each fit accurately by a power law; the indices of these power laws decrease 
with distance from the nucleus. B2 0326+39 also requires a discontinuity in emissivity 
to in order to fit a region with several bright knots of emission. In B2 1553+24, the 
sudden brightening of the jet can be explained by a combination of rapid expansion of 
the jet and a continuous variation of emissivity. The magnetic fields in both objects 
are dominated by the longitudinal component in the high-velocity regions close to 
the nucleus and by the toroidal component further out, but B2 0326+39 also has a 
significant radial component at large distances, whereas B2 1553+24 does not. Simple 
adiabatic models fail to fit the emissivity variations in the regions of high velocity but 
provide good descriptions of the emissivity after the jets have decelerated. Given the 
small angle to the line of sight inferred for B2 1553+24, there should be a significant 
population of similar sources at less extreme orientations. Such objects should have 
long (>200kpc), straight, faint jets and we show that their true sizes are likely to have 
been underestimated in existing images. 

Key words: galaxies: jets - radio continuum: galaxies - magnetic fields - polarization 
- MHD 



1 INTRODUCTION 

Evidence that the jets in low-luminosity, FRI 
iFanaroff fc Rilevl Il974l) radio galaxies are initially rel- 
ativistic and decelerate on kpc scales has mounted in 
recent years. FRI sources are thought to be the side-on 
counterparts of BL Lac objects, in which relativistic motion 
on p arsec scales is well-established llUrrv fc Padovanil 
Il995l) . Super luminal motions have been seen on milliarc- 
second scales in several FRI jets (|Gtov^nmni^^L|j20Wj) 
and on arcsecond scales in M 87 jBhe^a"^n^u^Towerj| 
Il995l) . In FRI sources, the lobe containing the main 



* E-mail: jrc@astro.ox.ac.uk 



(brig hter) jet is less de polarized than the counter-jet lobe 
llMorganti et alJ Il997bl) . This can be explained if the 
main jet points toward the observer, suggesting th at the 
brigh tness asymmetry is caused by Doppler beaming llLaind 
1988). T he asymmetry de creases with distance from the 
nucleus llLaing et al.l Il999l) . implying that the jets must 
decelerate. Self-consistent models of the deceleration of 
relativistic jets by injection of matter lost from stars or 
entrained from the surrounding galactic atmosphere have 
been calculated bvlBicknelll Il994) . iKomissarovl Jl994l) and 
iBowman. Leahy fc Komissarovl il99rJ) . 

lLaing fc Bridle! j2002al hereafter LB) fit VLA images 
of the radio jets in the nearby FRI radio galaxy 3C 31 
using a sophisticated model to reproduce the observed 
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total and polarized emission within 30 arcsec of the nu- 
cleus, where the jets are straight. They parameterized the 
three-dimensional distributions of velocity, emissivity and 
magnetic-field structure, calculated the brightness at each 
point within the jets in Stokes /, Q and U, accounting for 
the effects of relativistic aberration, and integrated along 
the line of sight to reproduce the expected distributions on 
the sky. They concluded that the jets in 3C 31 could be ac- 
curately modelled as intrinsically identical, axisymmetric, 
antiparallel, decelerating, relativistic flows, with locally ran- 
dom but anisotropic magnetic fields. Optimization of the 
model parameters placed tight constraints on the geome- 
try, vel ocity, emissivity and field structure. lLaing fcB ridlc 
(2002b) used this velocity field, together with estimates of 
the external pressure an d density from Chandra observations 
jHardcastle et all200 j) in a conservation-law analysis based 
on that of Bickncll ( 1994). They showed that there are self- 
consistent solutions for jet deceleration by injection of ther- 
mal matter and derived the variations of pressure, density, 
Mach number and ent rainment rate along the jets. Finally, 
lLaing fc Bridld i2004t) developed models of adiabatic, rela- 
tivistic jets with velocity shear and applied them to 3C31. 
They demonstrated that such models provide a reasonable 
description of the emissivity and magnetic-field variations 
at large distances from the nucleus but fail closer in, and 
inferred that significant reacceleration of relativistic parti- 
cles is requi red precisely where X-r ay synchrotron emission 
is observed iHardcastle et af]|2002l) . 

In the present paper we apply a modified version of LB's 
model to the jets of two FR I radio galaxies: B2 0326+39 
and B2 1553+24. Our principal aim is identical to that of 
LB: to estimate the distributions of velocity, emissivity and 
magnetic-field structure without introducing preconceptions 
about the (poorly understood) internal physics. By studying 
different objects we hope to be able to improve the range 
of physical scales we are able to probe, to identify which 
intrinsic features are common to all FRI jets and which 
vary from object to object and to assess the dependence of 
the jet structure on power and environmental conditions. 

Section [5] presents our new, deep, full-synthesis VLA 
images. In Section |3J we outline the model, emphasizing 
the parameterizations of the geometry, velocity, emissivity 
and magnetic field which differ from those used by LB. Sec- 
tion 0] compares our best-fitting models with the observed 
data in a variety of ways to show the features that we are 
able to reproduce as well at those we cannot. In Section^] 
we present the velocity, emissivity and field structures of 
the best-fitting models. Section HJ investigates whether the 
magnetic-field structure and emissivity are consistent with 
the assumptions of flux freezing and pure adiabatic energy 
loss and examines the idea that the jets reaccelerate. We 
then consider the appearance of the model for B2 1553+24 
at large angles to the line of sight and the implications for 
the detectability of side-on counterparts. Finally, we briefly 
compare the models for 3C 31, B2 0326+39 and B2 1553+24. 
Section [7| summarizes our conclusions and outlines possible 
avenues for further work. 

We adopt a Hubble constant, Ho = 70kms~ 1 Mpc _1 
throughout and define spectral index a in the sense S(v) oc 
v~ a . We use the notation p = (Q 2 + U 2 ) 1/2 /I for the degree 
of linear polarization. 



2 OBSERVATIONS AND IMAGES 

2.1 Object selection 

Our models assume that the jets are intrinsically identical, 
antiparallel, axisymmetric stationary flows, so any apparent 
differences must be due entirely to relativistic effects. Some 
important selection criteria are therefore imposed on our 
choice of source: 

(i) Only straight parts of the jets can be modelled and it 
is therefore important that the first bends in both the main 
and counter-jets be as far from the nucleus as possible. 

(ii) The jet and counter-jet must leave the nucleus in 
strictly antiparallel directions. 

(iii) We derive the sidedness ratio by dividing a total- 
intensity image of the jets by a copy of itself rotated through 
180° about the nucleus. The ratio is in the sense main 
jet/counter-jet and cannot be <1 over a significant area if 
the asymmetries are to be attributed entirely to relativis- 
tic effects (we expect some local fluctuations due to small- 
scale structure). This excludes a few objects, in particular 
those where t he counter-jet is w ider than the main jet (e.g. 
B2 0755+379; iBondi et al.ll2000D . 

(iv) There cannot be any other radio features (e.g. lobes) 
that could be confused with the emission from the jet. 

(v) The jets must be bright enough to produce a signifi- 
cant polarized signal from the region to be modelled. This is 
required in order to break the degeneracy between velocity 
and angle to the line of sight which is inherent in the total 
intensity (Section l3,21 . 

We selected two suitable targets from the well studied 
sample of nearby radio galaxies identified with sourc es in 
the B2 catalogue iColla et alJll975t iParma et al]ll987f) . 

2.2 B2 0326+39 

B2 0326+39 is a typical FRI radio sou rce associated wit h 
a bright elliptical galaxy at z = 0.0243 jMiller et alJl2002>) . 
Its luminosity at 1.4 GHz is P 1A = 2.2 x 10 24 WH z" 1 . 5- 
GHz images from the WSRT JParmall982ft and VLA feridld 
1982) first resolved the two fairly symmetrical, antiparal- 
lel jets. Comprehensive observations at multiple frequencies 
and resoluti o ns be tween 26 and 0.6 arcsec were made by 
iBridle et alJ 1^)91), who showed that the main (western) 
jet is much brighter than the counter-jet for the first 4 arc- 
sec from the nucleus and that the change from one-sided 
to two-sided structure is accompanied by a flip in the po- 
larization from parallel to perpendicular apparent magnetic 
field. The surrounding lobe emission is resolved out in im- 
ages with FWHM 2 arcsec. The mean spectral index of 
the jets close to the nucleu s is a = 0.55 between 1.4 and 

5.0 GHz iBridle et al.lll99ll) . 

ROSAT PSPC observations (ICanosa et al.ll999ft reveal 
a point source with a flux of 36 njy at 1 keV and an extended 
atmosphere well fitted by a beta model with form factor 
Patm = 0.35 and core radius r c = 30kpc. 

2.3 B2 1553+24 

B2 15 53+24 is a more distant (z = 0.0426; IColla et all 
Il975l) and less luminous (P XA = 4.7 x 10 23 WHz _1 ) 
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FRI radio galaxy. V LA observations at 1.4 and 4.8 GHz 
iMorganti et al Il98fl revealed two straight, well-collimated 
jets aligned N W - SE with very little surrounding lobe emis - 
sion (see also lStocke fc Burnslll98'7l Ide Ruiter et ai~lll993t) . 
Faint extended emission asso ciated with the brig hter jet is 
visible on the NVSS, however dCondon et al.119981) . The jets 
are similar at large distances from the bright compact nu- 
cleus but become very asymmetric in the central 5 arcsec 
Their mean spectral index is a = 0.60 jMorcanti^^al 
H587 ) . The bright radio, optical and X-ray co re lLaiTt^tHd 



Table 1 . Record of VLA observations 



1999HCapetti et alJ200otlCanosa et al.ll999h . the very large 



side-to-s ide asymmetries between the two jets close to the 
nucleus llLaing et al Il999h and the detection of optical emis- 
sion from the base of the main jet using HST iParma et al.l 
2003) all indicate that the jets are close to the line of sight. 



2.4 Observations and data reduction 

The observations were made using the VLA with a centre 
frequency of 8.46 GHz and a 100-MHz bandwidth in all four 
standard configurations for B2 0326+39 and in A, B and 
C configurations for B2 1553+24 (Table EJ. Either 3C48 or 
3C 286 was used as the primary calibrator for amplitude and 
one of 3C 138 or 3C 286 was observed to set the E- vector po- 
sition angle. The data for each configuration were reduced 
separately in the AIPS software package using standard 
techniques of calibration, imaging and self-calibration. The 
datasets were then concatenated using the AIPS task DB- 
CON starting with the widest configurations (A and B), first 
adjusting the core flux of the smaller-configuration dataset 
to match that observed with the larger configuration and im- 
aged at matched resolution. A further iteration of phase self- 
calibration was done after each combination. This procedure 
accounts for any core flux-density variation between observa- 
tions. CLEANed images were produced at three resolutions: 
full (FWHM 0.25 arcsec) and tapered to two lower resolu- 
tions giving better signal-to-noise ratio in the outer parts 
of the jets. I images were also produced using maximum- 
entropy deconvolution, first subtracting a point source at the 
position of the core. After the deconvolution step, the point 
source was added back in and the image convolved with the 
same truncated Gaussian beam as used in the CLEANed 
images. The maximum-entropy and CLEANed images were 
consistent apart from some low-level artefacts clearly associ- 
ated with the latter. We therefore show only the maximum- 
entropy I images. The Q and U images were CLEANed. 
Noise levels for / and Q/U images are given in Tableland 
are close to theoretical limits. All of the displayed polar- 
ization vectors have been rotated by 90° from the direction 
of the E- vector at 8.4 GHz and therefore represent apparent 
magnetic field. No high-resolution, multifrequency studies of 
Faraday rotation have yet been made for the two sources, so 
we have not corrected for this effect. The esti mates of rota- 
tion m easur e from two-frequency data giv en bv lBridle et al.l 
Jl99lf) and iMorganti et al.1 <1987t Il997al) and direct com- 
parison of their higher-frequency images with ours implies 
that the magnetic field directions inferred from our data are 
in error by <2° due to Faraday rotatio n. No corrections for 
Ricean bias JWardle fc Kronberdll974h have been made in 
deriving the images of degree of polarization, p, presented in 
this paper, but they have all been blanked at low signal-to- 
noise ratio as indicated in the figure captions, so the bias is 



Object 


Configuration 


Date 


Integration 
time (min) 


B2 0326+39 


A 


1999 Sep 18 


604 




B 


1999 Dec 19 


604 




C 


2000 May 14 


131 




D 


2000 Aug 17 


41 


B2 1553+24 


A 


1999 Sep 18 


603 




B 


1999 Dec 19 


321 




C 


2000 May 14 


107 



Table 2. Image resolutions and noise levels, crj is the off-source 
noise level on the I image; crp the average of the noise levels for 
Q and U. 



Object FWHM rms noise level 

(arcsec) (/^Jy / beam area) 

<7J <7p 



B2 0326+39 


1.50 


7.6 


6.2 




0.50 


5.7 


6.0 




0.25 


5.8 


6.S 


B2 1553+24 


1.50 


5.2 


5.8 




0.75 


5.7 


5.7 




0.25 


6.0 


6.5 



negligible. In any case, the fitting procedure described below 
(which uses the Q and U images directly) is unaffected. 



2.5 Source description 

The positions and 8. 4- GHz flux densities of the cores found 
from our A-configuration images are given in Table |3 



2.5.1 B2 0326+39 

Figs da) - (c) show the jets of B2 0326+39 at three differ- 
ent resolutions. In panel (a), the jets can be seen over their 
full extent. They are initially straight and antiparallel, al- 
though both bend slightly (by ~3°) about 25 arcsec from 
the nucleus. At a distance of ~22 arcsec the ridge line of the 
main jet starts to oscillate slightly. Panel (b) shows the in- 
ner 45 arcsec of both jets at 0.5-arcsec resolution: the jets 
appear very different from each other within ~ 10 arcsec of 
the nucleus but become almost indistinguishable further out. 
This is also illustrated in Fig^d), which shows the ratio of 

Table 3. Core flux densities and positions 



Object 


Core flux 


Core Position (J2000) 




(mJy) 


RA 


Dec 


B2 0326+39 


71.4 


03 h 29 m 23.887 3 


39°47'31.96" 


B2 1553+24 


63.5 


15 h 56 m 03.912 s 


24°26'52.92" 
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Figure 1. Montage showing the large-scale structure and jets of B2 0326+39 at 8.4 GHz. Panels (a) - (c) have North at the top; in the 
remaining panels, the images have been rotated slightly so that the jet axis is horizontal, (a) 5.7arcmin (170 kpc) East- West field at 
1.5 arcsec (0.75 kpc) resolution, (b) 1.5arcmin (45 kpc) East- West field at 0.50 arcsec (0.25 kpc) resolution, (c) 38 arcsec (19 kpc) East- 
West field at 0.25 arcsec (0.13 kpc) resolution, (d) The jet/counter-jet intensity ratio (defined as in Section 12.11 for the straight region of 
the jet within 22 arcsec of the nucleus at 0.50 arcsec (0.25 kpc) resolution. The grey-scale range (0 - 5) is indicated by the labelled bar. (e) 
Grey-scale of total intensity with superimposed polarization vectors for the straight region of the jets within 22 arcsec of the nucleus at 
0.50 arcsec (0.25 kpc) resolution. The lengths of the vectors are proportional to the degree of polarization p and their directions are those 
of the apparent magnetic field, (f) As in panel (e), but showing the inner 6 arcsec of the western jet at 0.25 arcsec (0.13 kpc) resolution. 
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1 D ancsec fS.B kpc) 



Figure 2. Montage showing the jets of B2 1553+24 at 8.4 GHz. Panels (a) - (c) show grey-scales of total intensity and have North at 
the top; the remaining images have been rotated so that the jet axis is horizontal, (a) lllarcsec (98 kpc) East- West field at 1.5arcsec 
(1.3 kpc) resolution, (b) 53arcscc (47 kpc) East- West field at 0.75arcscc (0.66 kpc) resolution, (c) 22.5arcsec (20 kpc) East- West field 
at 0.25arcsec (0.22 kpc) resolution, (d) The jet /counter- jet intensity ratio for the inner (straight) section within 30arcscc (26 kpc) of 
the nucleus at 0.75arcsec (0.66 kpc) resolution. The labelled wedge indicates the grey-scale range, (e) A grey-scale of the total intensity 
with superposed polarization vectors, showing 30arcsec from the nucleus of each jet at 0.75arcsec (0.66 kpc) resolution. The superposed 
vectors indicate the degree of polarization p (on a scale given by the labelled bar) and the direction of the apparent magnetic field, (f) 
As in panel (e), but for the inner 5arcsec of the main jet at 0.25arcsec (0.22 kpc) resolution. 
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the main and counter-jet intensities at the same resolution. 
Panel (c) shows the jets at our highest resolution, 0.25 arcsec. 
The brightest region of the main jet is resolved into three 
knots embedded in more extended emission. The nucleus has 
an extension pointing in the direction of the main jet whose 
brightness drops rapidly with distance from the nucleus; it 
can only be seen out to «1 arcsec. No emission is visible at 
0.25-arcsec resolution between this point and the first bright 
knot of emission ~2 arcsec from the nucleus. 

The polarization structure of B2 0326+39 is shown at 
resolutions of 0.50 and 0.25 arcsec in Figs UTe) and (f), 
respe ctively. As is usual in FRI objects feridle fc Perlevl 
1984), the main jet close to the nucleus is polarized with an 
apparent field parallel to its axis. Further out in the main jet 
and in all parts of counter-jet with sufficient signal-to-noise 
ratio, the apparent field is perpendicular to the axis. In this 
object, unlike 3C31 and B2 1553+24, no parallel field edge 
is seen. 

We model the jets within 22 arcsec of the nucleus, where 
their ridge lines are straight. 

2.5.2 B2 1553+24 

Figs |3 a) - (c) show the jets of B2 1553+24 at three differ- 
ent resolutions: 1.5, 0.75 and 0.25 arcsec. The jets are very 
straight with only a slight bend in the main jet about 45 arc- 
sec from the nucleus. Both jets gets steadily fainter further 
out, and are lost in the noise approximately 65 and 45 arcsec 
from the nucleus for the main and counter-jets respectively 
in our lowest (1.5 arcsec) resolution image. The counter-jet is 
invisible within 3 arcsec of the nucleus; in contrast, the main 
jet is very bright there. Further out, the counter-jet bright- 
ens, reaching a maximum at 7 arcsec from the nucleus. From 
this point outwards, the main and counter-jets appear essen- 
tial identical, but with the main jet approximately twice as 
bright. This is also illustrated in the image of jet/counter- 
jet intensity ratio (Fig.|2]i). Initially (within 1.5 arcsec) the 
main jet is well collimated, as are both jets at distances 
^7 arcsec. Both jets flare dramatically between these regions 
of good collimation. 

The polarization structure of the jets in B2 1553+24 
is shown in Figs |5fe) and (f) at resolutions of 0.75 and 
0.25 arcsec respectively. The bright 4 arcsec of the main jet 
closest to the nucleus has an apparent magnetic field aligned 
along the jet axis. Further out, the field becomes perpendic- 
ular to the axis close to the centre of the jet and parallel to 
the outer isophotes near its edges. 

We model the jets out to 30 arcsec from the nucleus. 
This limit is set not by changes in jet direction, as in 3C 31 
and B2 0326+39, but by signal-to-noise, particularly in po- 
larization. 



spectrum n{E)dE = no-E - di? (corresponding to a fre- 
quency spectral index a) with an isotropic pitch-angle distri- 
bution. The maximum degree of linear polarization is then 
po = (3a + 3)/(3a + 5). We assume that the 1.4 - 5 GHz 
spectral indices given in Sections 12.5.11 and 12.5.21 also apply 
at 8.4 GHz. 

(iii) The magnetic field is tangled on small scales, but 
anisotropic (the reas ons for taking the field to be of this 
form are discussed bv lLain dll98ll . lBegelman et alJll984l and 
LB). 

We define f3 = v/c, where v is the flow velocity. V — (1 — 
/3 2 ) -1 / 2 is the bulk Lorentz factor. 

3.2 Overview 

We aim to determine the velocity and angle to the line 
of sight independently by comparison of emission from the 
main and counter-jets and the modelling of linear polariza- 
tion is essential to our technique. In order to illustrate this 
point, we consider the simple example of cylindrical, antipar- 
allel jets with constant velocity /3c at an angle 9 to the line 
of sight. If we consider only total intensity, the ratio of the 
flux densities per unit length for the main and counter-jets, 
7j/7 c j, does not uniquely determine the velocity: 

Ij_ / l + /3cos6> \ 2+ ° 

Jcj ~~ \l~f3cos9J 

for emission which is isotropic in the frame of the jet flow. 
In order to break the degeneracy between j3 and 9, we use 
the linear polarization. The relation between the angles to 
the line of sight in the rest frame of the flow, 9' and in the 
observed frame, 6, is: 

sin^' = [r(l — /3cos #)] -1 sin 6 (main jet) 
sin#cj = [r(l + /3cos #)] _1 sin 9 (counter-jet) 

The observed polarization is in general a function of 0' if the 
field is anisotropic. If we know the field structure a priori, 
then we can solve explicitly for j3 and 9. We take the example 
of a field which is disordered on small scales but confined to 
a plane perpendicular to the jet, with equal rms along any 
direction in the plane. In this case, there is no variation of 
the degree or direction of polarization across the jet. For 
a = 1, the total and polarized flux densities per unit length 
in the emitted frame are: 

/' = K(l + cos 2 $') 
P' = poK sin 2 0' 



3 THE MODEL 
3.1 Assumptions 

Our fundamental assumptions are those of LB: 

(i) The jets may be modelled as antiparallel, axisymmet- 
ric, stationary, laminar flows. 

(ii) They contain relativistic particles with an energy 



where K is a constant (lLaindll980l.ll98 J) . The ratios of ob- 
served total and polarized intensity for the jet and counter- 
jet are: 

3 



[r(i 



/3cos6>)]" 2 sin 2 6> 



[r(l + /?cos6>)]- 2 sin 2 
1 + /3cos# x " 

1 -/?COS0 
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These equations c an be solved numerically for f3 and 8 (cf. 
iBondi et al.ll200ol) . [The specific field configuration is a rea- 
sonable approximation for the outer regions of the jets in 
B2 0326+39 (Section EH ]. 

In general, we must fit the field configuration and there- 
fore need to introduce two additional parameters to describe 
its anisotropy (Section 13,51 . These also determine the vari- 
ation of polarization transverse to the jet axis, and can be 
estimated independently of the velocity and angle if the jets 
are well resolved in this direction. 

Following LB, we model the jets using simple parame- 
terized forms for the jet geometry, velocity, emissivity and 
magnetic-field structure. We find that the functional forms 
used by LB to model variations along the jets are unable 
to fit the new data. In order to produce models that fit the 
data accurately, we are forced to decouple the parameter- 
izations of velocity, emissivity and magnetic field from the 
three regions defined by the geometry in LB. This provides 
much more freedom in our models, but also allows simpler 
functional forms to be used. 

LB examined two forms for transverse variations of ve- 
locity and emissivity: spine /shear-layer (SSL) and Gaussian. 
The former gave a slightly better fit to the data for 3C31, 
at the cost of introducing extra parameters. We find that 
the off-axis variations of velocity and emissivity are less well 
constrained in our objects and that the difference between 
SSL and Gaussian models is insignificant. In what follows, 
we use only the simpler Gaussian models. 

3.3 Geometry &: Coordinate Systems 

The geometry of the model jet is sketched in Fig. [3] It is 
divided into two parts: a flaring region where the jet first 
expands rapidly and then recollimates and an outer region 
of constant opening angle. The mathematical description of 
the edge of the jet is as given by LB, but with two important 
changes. First, the conical outer region is not centred on the 
nucleus but at a distance A behind the nucleus along the jet 
axis, allowing the jet to recollimate much more than is seen 
in 3C31 (Fig. |3J . This modification is required to describe 
many FRI jets, including both B2 0326+39 and B2 1553+24. 
Second, we find no evidence for separate conical inner re- 
gions in the B2 sources (the jets are very narrow and poorly 
resolved close to the nuclei, so their collimation behaviour 
is poorly constrained there) . We therefore extend the flaring 
region to the base of the jet. An inner region analogous to 
that used by LB to model 3C 31 is defined purely by the 
emissivity profile (Section l3.6|I .The distance of the outer- 
most isophote from the jet axis in the flaring region is fit by 
a cubic polynomial, as in LB, but both it and its derivative 
are constrained to be zero at the nucleus (£i = and n = 
in the notation of LB). 

Four parameters define the geometry of the model. Of 
these, the angle to the line of sight, 8, the opening angle of 
the conical outer region, £o, and the distance along the axis 
from the nucleus to the boundary between flaring and outer 
regions, ro, are defined essentially as in LB. Because the 
vertex of the outer region is no longer at the nucleus, how- 
ever, the spherical boundary surface is centred on the vertex 
and has radius ro + A (this reduces to the case considered 
by LB if A = 0). The fourth defining parameter is xo, the 
radius of the jet edge at the outer boundary (Fig. [3J . We 



choose to optimize xo rather than A to avoid cross-coupling 
between variables: xo is determined by the position of the 
outer isophote at the boundary and £o by the expansion of 
the jet in the outer region, whereas A affects both quantities. 
The variables are related by (ro + A) sin£o = xo- 

The velocity, emissivity and magnetic-field parameteri- 
zations defined below are expressed in a streamline coordi- 
nate system (p, s) where the streamline index s is constant 
for a given streamline and p increases monotonically with 
distance along it. The distance of a streamline from the jet 
axis is: 

x(z,s) = a2{s)z 2 + as(s)z 3 (flaring region) 
x(z,s) = (z + A) tan(£os) (outer region) 

where z is the distance from the nucleus along the axis. In 
the outer region, s = f/fo, where £ is the angle between 
the flow vector and the jet axis. For the flaring region, 02 (s) 
and a3(s) are defined by the conditions that x(z,s) and its 
derivative with respect to z, x'(z,s), are continuous at the 
boundary between the two regions, s therefore varies from 
on-axis to 1 at the jet edge; the s — 0.5 streamline is plotted 
in Fig.|3] The form of x(z, s) in the flaring region is that used 
by LB with the simplifying condition x(0, s) = x'(0, s) = 
0. Parameterization of physical quantities in terms of the 
distance from the nucleus along a streamline creates undue 
complexity because of the cubic form of x(z, s). Instead, we 
use a coordinate: 

zro 

P (r + A)cos^ s)-A P<r ° 

Z + A A 

cos(£os) 

p increases monotonically along a streamline from at the 
nucleus and each streamline crosses the boundary into the 
outer region at p = ro. The coordinate system used by LB 
is the special case with A — 0. 

3.4 Velocity field 

The variations in sidedness radio along the jets in both 
sources (Figs^i,|5ji) suggest that the initial velocity needs 
to be high in order to account for the asymmetry as a con- 
sequence of Doppler boosting. Further from the nucleus, 
there is a transition to a regime where the sidedness ra- 
tios are low and approximately uniform. In B2 1553+24, the 
ratio remains constant, but significantly >1 for all distances 
>7arcsec from the nucleus, indicating a mildly relativistic 
but constant velocity. Two changes to the form of the veloc- 
ity profile used by LB to describe 3C 31 are required: 

(i) In the flaring region, the function used by LB describes 
an initial fast section and abrupt deceleration, just as needed 
for the B2 sources, but the deceleration is forced to occur 
immediately before the jet recollimates. This does not ap- 
pear to be general, so we need a functional form which is 
qualitatively similar to that used by LB, but which allows 
the location and length of the rapid deceleration region to 
be varied. 

(ii) The jet/counter-jet sidedness ratio in the outer region 
of 3C 31 continues to decline slowly with distance from the 
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Figure 3. Geometry of the model. The jet edge (s = 1 streamline) is indicated by the solid line and the s = 0.5 streamline by the dashed 
line. The arc shows the spherical boundary between the flaring and outer regions, centred at the vertex of the outer region. The three 
quantities used to define the jet geometry: ro, the distance of the outer boundary from the nucleus; xq, the width of the jet at the outer 
boundary and £o> the opening angle of conical outer region are marked. 



nucleus, so LB used an exponential form for the velocity 
profile there. This will not fit the large regions of constant 
sidedness ratio in the B2 sources, especially in B2 1553+24. 

We therefore divide the on-axis velocity profile into three 
regions: (a) approximately constant, with a high velocity 
close to the nucleus; (b) a linear fall-off and (c) roughly 
constant, but with a low velocity at large distances. Example 
profiles are given later (Figs UTTb and 118b ). The profile is 
defined by four free parameters: the distances of the two 
boundaries separating the three regions, pp y and p@ , and 
the characteristic inner and outer velocities /3i and /3o. 

The velocity along any off-axis streamline is calculated 
using the same expressions but with inner and outer veloci- 
ties f3i exp(— s 2 Invi) and f3o exp(— s 2 ln«o), respectively, i.e. 
with a truncated Gaussian transverse profile falling to frac- 
tional velocities Vi and vo at the edge of the jet in the inner 
and outer regions, respectively. This differs slightly from the 
separable function used by LB, but gives essentially identical 
transverse profiles. 

The full functional forms for the velocity field /3(p,s) 
axe given in Table 0] [Note that the exponential terms are 
introduced purely to maintain continuity in the acceleration 
near the boundaries and their precise form has no significant 
effect on the brightness distribution]. The constants ci, C2, 
C3 and C4 are defined by the values of the free parameters 
and the conditions that the velocity and acceleration are 
continuous at the two boundaries. 

3.5 Magnetic Field Structure 

We define the rms components of the magnetic field 
to be (Bf) 1 /' 2 (longitudinal, parallel to a streamline), 
{B 2 ) 1 ^ 2 (radial, orthogonal to the streamline and outwards 
from the jet axis) and (B 2 ) 1 ^ 2 (toroidal, orthogonal to the 
streamline in an azimuthal direction). The magnetic-field 
structure is parameterized, as in LB, by the ratio of rms ra- 
dial/toroidal field, j(p,s) = (B 2 y^/{B 2 ) 1/2 and the longi- 
tudinal/toroidal ratio k(p,s) = {Bf ) 1/2 /{B 2 } 1/2 . The main 
difference between our model and that of LB is again that 
the characteristic distances are defined independently from 
those for other quantities. Two fiducial distances, pb 1 and 
Pb , are used. For p < pb 1 and p > pB , the field ratios have 



constant values, with linear interpolation between them for 
pBi *S p ^ pB - We find no evidence for any transverse 
structure in either of the field ratios. In 3C 31, there are sig- 
nificant variations in the radial/toroidal field ratio between 
the centre and edge of the jet, and the model of LB is there- 
fore more complex than that used here. The functional forms 
assumed for the field ratios are again given in Table 2] 

3.6 Emissivity 

We write the proper emissivity as e(p,s)f(p,s), where e is 
the emissivity in / for a magnetic field B = {B? + B 2 + 
B 2 ) 1 / 2 perpendicular to the line of sight. / depends on field 
geometry: for /, ^ / ^ 1 and for Q and U —po ^ / ^ +po- 
Following LB, we refer to e, loosely, as 'the emissivity'. For 
a given spectral index, it is a function only of the rms total 
magnetic field and the normalizing constant of the particle 
energy distribution, e oc noB +OL . 

The description of the emissivity is similar to that given 
by LB: the on-axis profile is divided into distinct regions, 
each with a power-law profile. As for the velocity, however, 
the boundaries between regions are defined without refer- 
ence to the geometry. In order to differentiate between sig- 
nificant changes in the emissivity profile and small-scale, 
stochastic features (which will occur at different places in 
the main and counter-jets), we require evidence that the 
same emissivity structure is present in both jets and that 
it extends over many beam areas. Five regions are required 
(although only four are used for each source) compared with 
the three used by LB. The reasons for the increase in the 
number of regions are: 

(i) In B2 0326+39 the bright region of the main jet be- 
tween 2 and 5 arcsec from the nucleus (Fig. \^>) is composed 
of three knots of emission of approximately equal brightness 
with a sharp falloff at each end; there is a much fainter equiv- 
alent in the counter-jet. Two regions are therefore required 
rather than the one used in LB: a relatively flat section aver- 
aging the intensity of the knots and a steep section modelling 
the rapid falloff at the end of the last knot. The rapid in- 
crease in brightness at the innermost knot is modelled using 
a discontinuity in the emissivity profile as in LB. 

(ii) A further region is required to model the outer parts 
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Table 4. Functional forms of the velocity 0, emissivity e, radial/toroidal and longitudinal/toroidal magnetic-field ratios j and A; in the 
streamline coordinate system (p, s). Column 4 lists the parameters which may be optimized, for comparison with Table© 



Quantity Functional form Range Free parameters 



pi cxp( — s In v\ ) — /3q exp( — s ln^o) 
10 



Velocity field 

exp[ci(p — pvi)] P < Pv! Distances p vi , p vo 

C2 + C3p pvi ^ p < pv Velocities /3i, /3o 

/3i — ® 1 cxp ^ — ° lnt,1 ^ 10 ^° cxp ^ — s ln "°) exp[c4(pv — />)] p > Pv Fractional edge velocities ui, jjq 



Emissivity 



d 2 



cxp 



: In 



-B 3 



^ J exp(-s 2 lne ) 



exp(— s 2 Ineo) 



) exp(-s 2 lne ) 



p ^ p ei Distances p ci ,p C2 , Pe 3 ,Pe 4 



Pci < p $5 p C2 Indices Si, i?2, E3, E4, E5 



pc 2 < p pc 3 Fractional increase g 



Po 3 < P < Pe 



P > Pe 4 



i(p,s) Jl 

JO 



Radial/toroidal field ratio 

P ^ PBi Distances p Bl , Pb 

PBi < P < PB Ratios ji, jo 
P ^ Pb 



k(p,s) ki 



* + (*»-*0(£^) 



fen 



Longitudinal/toroidal field ratio 

P ^ PBi Ratios fci, fco 

PBi < P < PB 

P > PB 



of both jets in B2 1553+24 (we will show that this is on a 
much larger physical scale than the entire modelled regions 
of 3C 31 and B2 0326+39; it is therefore reasonable that it 
shows different emissivity behaviour). 



The emissivity is continuous across all of the boundaries 
except the innermost (p = p ci ), where a jump by a factor of 
g is allowed. The part of the jet with p < p ci is analogous to 



the inner region defined by geometry, emissivity and velocity 
for 3C31 (LB). 

Off-axis, the profile is multiplied by a factor 
exp[— s 2 lne(p)], as in LB, so that e(p) is the fractional value 
of the emissivity at the jet edge. e(p) has a constant value 
eo for p > p 02 and varies linearly through region 2 from e\ 
at p ci to eo at p C2 . In the inner region the jet is too narrow 
to constrain any transverse profile and we set e(p) = 1. 
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Table 5. The number of independent points and the reduced \ 2 
values for our final optimized models. 



Object 


Independent 


Reduced \ 2 




points 




B2 0326+39 


3614 


1.32 


B2 1553+24 


3722 


1.13 



The full description of the emissivity distribution e(p, s) 
is given in Table 0] 



3.7 Modelling Procedure 

The modelling procedure is that described by LB, with one 
small refinement. In order for the numerical integrations to 
converge, it is necessary to determine appropriate limits, 
which are where the line of sight intersects the jet surface. 
If 9 is small, a line of sight may cross the jet boundary more 
than twice, for example in passing through both the inner 
part of the jet and the shoulder where the jet recollimates. 
This possibility is required for B2 1553+24 and is now in- 
cluded. 



3.8 Fitting and optimization 

We use x 2 as a measure the goodness of fit. The noise level 
for each resolution and Stokes parameter is estimated as 
l/y/2 times the rms difference between the image and a copy 
of itself reflected across the jet axis, as in LB. This quantity 
is always significantly larger than the off-source noise level. 
It is a measure of the local deviation from mirror symmetry, 
which limits our ability to fit the data using an axisym- 
metric model, and includes a contribution from any residual 
deconvolution errors. It is a lower limit to the correct noise 
level, as it excludes any mirror-symmetric component. We 
therefore expect \ 2 to be overestimated. 

X 2 is calculated using the high-resolution (0.25 arcsec) 
images where they have adequate signal-to-noise ratio (i.e. 
close to the nucleus) and the lower-resolution images else- 
where. x 2 values for /, Q and U are summed. In a slight 
modification to the algorithm used by LB, the sizes of the 
regions over which x 2 ls evaluated from the high-resolution 
data may be different for total intensity and linear polar- 
ization. For B2 0326+39 we fit to the 0.25-arcsec / within 
4.5 arcsec, further out, we use the 0.5-arcsec image; for the 
Q and U images only the 0.5-arcsec images are used. For 
B2 1553+24 we fit the inner 4 arcsec of the 0.25-arcsec I, Q 
and U images and the 0.75-arcsec images at larger distances. 
The sampling grid is chosen so that all estimates of x 2 are 
independent and the core is always excluded. 

An initial starting model is generated by adjusting a 
number of key parameters by hand to match recognizable 
features of the jet. Their values are then fixed and the re- 
maining parameters are optimized using the downhill sim- 
plex method of Neld er fc Mead implem ented using Numeri- 
cal Recipes routines JPress et al.lll992l) to minimize x • 



3.9 Uniqueness and errors in parameter 
estimation 

The uniqueness of such a complex model is a concern, espe- 
cially given that the downhill simplex algorithm is not guar- 
anteed to converge to a global x 2 minimum. In addition, the 
parameters are clearly not independent, and coupling be- 
tween them leads to some indeterminacy. We are confident 
of the reliability of our models, however, for the following 
reasons: 

(i) We fit to deep, well-resolved images in I, Q and U. 
These provide a large number of independent data points 
with good signal-to-noise ratio. 

(ii) We have experimented with a wide range of starting 
conditions and have been unable to find satisfactory solu- 
tions which differ significantly from those presented here: 
models even remotely resembling the data are hard to find. 

(iii) Whilst changes in some parameters can, in part, be 
offset by the variation of others, the extent to which this can 
occur is strictly limited by the fact that most parameters 
only affect the fit in limited regions and are well constrained 
by gross features in the brightness distributions. 

We derive rough uncertainties, as in LB, by varying individ- 
ual parameters until the increase in x 2 f° r total intensity or 
linear polarization in either the flaring or the outer region 
corresponds to the formal 99% confidence level for indepen- 
dent Gaussian errors. These estimates are crude (they ne- 
glect coupling between parameters), but in practice give a 
good representation of the range of qualitatively reasonable 
models. 



4 COMPARISON BETWEEN MODELS AND 
DATA 

4.1 x 2 values 

The values of reduced x 2 f° r the final optimized models of 
the two sources are listed in Tabled] Given that the "noise 
levels" are probably underestimated (Section l3.8H . these rep- 
resent very good fits. In the remainder of this section, we 
compare the models and data in detail. 

4.2 B2 0326+39 

4-2. 1 Total intensity images and profiles 

Figs 0| and [S] show model and observed contour plots and 
profiles along the jet axis for Stokes / at 0.50 and 0.25 arcsec 
resolution, respectively. The model and observed sidedness 
ratios are shown as grey-scale plots and on-axis profiles in 
Fig.© 

4-2.2 Fitted total intensity features 

Our optimized model successfully reproduces the following 
features of the total-intensity distribution in B2 0326+39, as 
illustrated in Figs ^ — |H] 

(i) The base of the main jet is initially bright but fades 
rapidly away from the nucleus and becomes invisible after 
1 arcsec. 
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Figure 4. A comparison of the model and data for B2 0326+39 in total intensity at 0.50 arcsec resolution for the inner 21 arcsec of 
each jet. (a) model contours; (b) observed contours. The contour levels in panels (a) and (b) are 1, 2, 4, 8, 16, 32, 64, 128, 256, 512 X 
20/iJy/beam area, (c) Profiles along the jet axis for the data (solid line) and model (dashed line). 



(ii) Apart from the bright base of the main jet, both jets (iv) A similar, but much fainter, peak in the brightness 

are very faint within 2 arcsec of the nucleus. of the counter-jet is also reproduced. 



(iii) The main jet has a very bright section between 2 
and 4.5 arcsec from the nucleus, brightening and fading very 
quickly at each end. 



(v) The main jet beyond 5 arcsec has a very flat profile, in 
contrast to the counter-jet which peaks at « 15 arcsec from 
the nucleus. 
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Figure 5. A comparison of the model and data in total intensity 
at 0.25 arcsec resolution for the inner 6 arcsec of the main jet 
in B2 0326+39. (a) model contours; (b) observed contours. The 
levels in panels (a) and (b) are 1, 2, 4, 8, 16, 32. 64, 128, 256, 
512x 25^tJy/beam area, (c) Profiles along the jet ridge-line for 
the data (solid line) and model (dashed line). 

(vi) The jets are initially well collimated, flare 5 arcsec 
from the nucleus and recollimate at «15 arcsec. 

(vii) The on-axis sidedness ratio falls over the first half of 
the modelled region to near unity at 11 arcsec. This value is 
maintained over the outer 10 arcsec. 
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Figure 6. A comparison of the model and observed jet/countcr- 
jet sidedness ratios at 0.50 arcsec for B2 0326+29. (a) Model grey- 
scale; (b) observed grey-scale. The grey-scale range for panels (a) 
and (b) is - 5, as indicated by the labelled bar. (c) Profiles along 
the jet ridge-line for the data (solid line) and model (dashed line). 



4-2.3 Polarization images and profiles 

Fig. |7]illustrates the observed and model polarization struc- 
tures of the jets in B2 0326+39 as grey-scales of the degree 
of polarization, p, vector plots representing p and the ap- 
parent magnetic field direction and profiles of p along the 
jet centre-line. There is little transverse variation of the de- 
gree or direction of polarization across the jets at large z 
(Figs [TJd and d) and the ridge-line profiles are quite noisy, 
so we also show profiles of p averaged across the jets from 
10 - 21.5 arcsec from the nucleus in Fig. |S] Vector plots for 
the inner regions of the main jet are shown on an expanded 
scale in Fig. fj The 0.25 arcsec resolution images have low 
polarized signal-to-noise ratio and are not shown. 

Individual transverse profiles of p are also noisy, so 
Fig. 1101 shows averages over the outer parts of the main and 



counter-jets (between 15 and 21 arcsec from the nucleus), 
where there is little longitudinal variation. 



4-2.4 Fitted polarization features 

Our model is successfully able to reproduce the following 
features in the polarization images (Figs I7 HTU1 : 

(i) The inner knot and the bright regions of the main 
jet within 5 arcsec of the nucleus have longitudinal apparent 
magnetic field. 

(ii) The bright part of the main jet is more polarized at 
its edges, forming a V-shaped structure in p (Figs an d 
b). 

(iii) In this region,the on-axis degree of polarization falls 
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Figure 7. A comparison of the model and observed polarization distributions at 0.50 arcsec resolution in the inner ±21 arcsec of the jets 
in B2 0326+39. (a) and (b): grey-scales of the degree of polarization, p. The grey level is indicated by the labelled bar at the top of the 
figure, (a) model; (b) data, (c) and (d): vector plots showing the degree and direction of linear polarization superimposed on grey-scales 
of total intensity. The lengths of the vectors are proportional to p (on a scale given by the labelled bar in panel c) and their directions 
are those of the apparent magnetic field, (c) model; (d) data, (e) Profiles of degree of polarization along the centre-line of the jet for the 
data (solid line) and model (dashed line). In all displays, the data and models have been blanked wherever the polarized power is < 3<Tp 
or the total intensity is < 5o7, using the values of off-source noise given in Table |21 
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Figure 8. Profiles of the degree of polarization averaged across 
the jets of B2 0326+39. Full line:data; dashed line: model. The 
plots cover the range 10 - 21.5 arcsec from the nucleus, where 
there is little variation across the jets (Fig. [TuV (a) counter-jet; 
(b) main jet. The resolution is 0.5 arcsec FWHM. 



B2 0326+39 




Figure 9. Vector plots of the degree and direction of linear po- 
larization superimposed on grey-scales of total intensity for the 
inner 6 arcsec of the main jet in B2 0326+39 at a resolution of 
0.50 arcsec. The vector lengths are proportional to the degree of 
polarization, p and their directions are those of the apparent mag- 
netic field, (a) model; (b) data. 



from 0.2 at the base of the main jet to 0.1 at 4 arcsec 

(Kg. no. 

(iv) The main jet has a region where no polarized signal 
is detected at the 3cr level between 5 and 9 arcsec, setting a 
limit of p < 0.15. The model distribution is consistent with 
this limit. In the counter-jet, no significant polarized signal 
is seen closer than 10 arcsec from the nucleus and none is 
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Figure 10. Transverse profiles of the degree of polarization 
in B2 0326+39 obtained by averaging over z between 15 and 
21 arcsec. Full line: data; dashed line: model. The observed profiles 
are derived by averaging over images of p made without blanking, 
in order to avoid biasing the results to high values. The resolution 
is 0.5 arcsec FWHM. 



predicted (the observational limits are not severe because of 
the low total intensity). 

(v) Beyond 10 arcsec from the nucleus, the jet is polarized 
such that the apparent magnetic field vector is transverse to 
the jet axis and no parallel-field edge is seen. 

(vi) Between 10 and 15 arcsec, the polarizations of both 
the jet and counter-jet increase monotonically before be- 
coming constant with p w 0.5 further out. [The apparent 
underestimation of p in the counter-jet at as 10 arcsec from 
the nucleus is in a region of low polarized signal and is prob- 
ably not significant]. 



4-2.5 Features that cannot be fitted well 

Our model cannot reproduce small-scale and/or non- 
axisymmetric features, the most prominent of which are the 
three knots of emission in the bright region of the main jet 
between 2 and 4.5 arcsec from the nucleus (Fig|jjj). The model 
correctly averages the brightness in this region into a smooth 
profile. In addition, there are a number of small, but signif- 
icant discrepancies between model and data: 

(i) The inner region of the counter-jet is observed to be 
brighter than that of the main jet (Fig^J)), an effect that 
is impossible to reproduce using intrinsically identical jets. 
This may also be caused by fluctuations on a small scale 
(e.g. a single knot in the counter-jet). 

(ii) The counter-jet first brightens significantly 1 arcsec 
further from the nucleus than the main jet and the minimum 
in the total intensity profile at the end of the flaring region is 
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Figure 11. A comparison of the model and observed total intensity for B2 1553+24 at 0.75 arcsec resolution. The plots show the inner 
±29 arcsec. (a) and (b): contour plots of total intensity with levels at 1, 2, 4, 8, 16, 32, 64, 128, 256, 512x 20/iJy/beam area, (a) model; 
(b) data, (c) Profiles along the jet centre-line for the data (solid line) and model (dashed line). 



2 arcsec further from the nucleus in the counter-jet (FigEt)- 
This is not reproducible by our model, for which the overall 
shape of the flaring region should be similar for the two jets, 
even if their brightnesses differ significantly. 

(iii) The counter-jet intensity decreases with z faster than 
the model predicts for z > 15 arcsec. This has the effect of 
increasing the sidedness ratio slightly in this region, an effect 
which is inconsistent with a monotonically decelerating flow 
(see Section lb.2l . 

(iv) The model predicts a degree of polarization which is 
slightly higher for the main jet than the counter-jet at large 
z. The converse is true, although the discrepancy is marginal 
(Figs|7J; andO§. 

(v) In the outer parts of the modelled region, p is pre- 
dicted to be higher on-axis than at the jet edges. This is 



marginally inconsistent with the data for the main jet, where 
p has little if any transverse variation (Fig. I1UI . 

4.3 B2 1553+24 

4-3. 1 Total intensity images and profiles 

Figs II II shows contours and longitudinal profiles of total in- 
tensity for the model and observed images at 0.75 arcsec 
resolution. Corresponding plots for the inner 5 arcsec of the 
main jet at 0.25-arcsec resolution are shown in Fig. 1121 (the 
counter-jet is invisible close to the nucleus on these images) . 
The model and observed sidedness ratios are compared in 
Fig. ESI 
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Figure 12. A comparison of the model and observed total inten- 
sity at 0.25 arcsec resolution for the inner 5 arcsec of the main jet 
in B2 1553+24. (a) and (b): contour plots with levels at 1, 2, 4, 
8, 16, 32, 64, 128, 256, 512 X 20A*Jy/beam area, (a) model; (b) 
data, (c) A profile along the jet centre-line showing the observed 
(solid line) and model (dashed line) values. 
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Figure 13. A comparison of the model and observed jet/counter- 
jet sidedness ratios at 0.75 arcsec for B2 1553+24. (a) and (b): 
grey-scale images of sidedness ratio, in the range - 5 as indi- 
cated by the labelled bar. (a) model; (b) data.(c) sidedness profiles 
along the centre-line of the jet showing data (full line) and model 
(dashed line). 



0.75-arcsec resolution image. After this point it brightens, 
reaching a peak ~7 arcsec from the nucleus. 

(v) At z « 6 arcsec in the main jet, the total intensity 
profile flattens before falling again from z « 8 arcsec. 

(vi) The jet/counter-jet sidedness ratio drops rapidly 
with z, reaching a constant value of ~2.2 at z fa 6 arcsec. 



4-3.2 Fitted total intensity features 

The following features of the jets in B2 1553+24 are fitted 
accurately by our optimized model: 

(i) The main jet initially brightens very rapidly away from 
the nucleus, peaking at 0.7 arcsec separation and then fading 
to give the appearance of a bright knot. 

(ii) This jet remains well-collimated for 1.5 arcsec before 
flaring abruptly. 

(iii) At z as 7 arcsec, both jets recollimate, maintaining 
an approximately constant width at large z. 

(iv) The counter-jet is not seen in the inner 3 arcsec of the 



4-3.3 Polarization images and profiles 

Fig |14l shows grey-scale images of the degree of polarization 
p and vector images showing the orientation of the appar- 
ent magnetic field at both 0.75 and 0.25 arcsec resolution, 
together with profiles of p along the jet axis. Transverse pro- 
files of the degree of polarization at two representative points 
in the jet and one in the counter-jet are shown in Fig 1151 



4-3-4 Fitted polarization features 

Our model is successfully able to reproduce the following 
features of the linear polarization in B2 1553+24; 
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Figure 14. A comparison of the model and observed polarization for B2 1553+24. Panels (a) — (e) show the inner ±29 arcsec at 0.75 arcsec 
resolution, (a) and (b): grey-scale images of the degree of polarization, p, in the range — 0.7, as indicated by the labelled wedge, (a) 
model; (b) data. Panels (c) and (d) are vector plots showing the degree of polarization and the apparent magnetic field direction. The 
vector lengths are proportional to p, on a scale indicated by the bar on panel (c). (c) model; (d) data, (e) Profiles of the degree of 
polarization along the centre-line of the jet for model (dashed) and data (full). The data and models have been blanked wherever the 
polarized signal is < 3<rp or the total intensity is < 5<tj, using the values of off-source noise as given in Table |21 Panels (f) - (j): as (a) - 
(e), but for the inner 5 arcsec of the main jet at a resolution of 0.25 arcsec. Note that the grey-scale range in panels (f) and (g) is — 0.5 
and therefore differs from that in panels (a) and (b). 



(i) Within 3.5 arcsec of the nucleus, the degree of polar- 
ization is approximately constant on-axis, with p m 0.15. It 
increases with distance from the axis, reaching ~0.3 at the 
edge of the jet. 

(ii) In this region the apparent field is orientated approx- 
imately parallel to the jet axis. 



(iii) At «4 arcsec from the nucleus the polarization on- 
axis in the main jet is very low. 

(iv) Between 4 and 8 arcsec from the nucleus the on-axis 
polarization in the main jet increases to p ~0.4; further out 
it remains approximately constant at that value. 
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Figure 15. Transverse profiles of the degree of polarization in 
B2 1553+24. (a) Profile at a distance of 1.7 arcsec from the nu- 
cleus in the main jet. The apparent field is parallel to the jet axis 
across the entire width, (b) and (c) Profiles at 7.5 arcsec from 
the nucleus for the main and counter-jets, respectively. In both 
jets, the apparent field is parallel to the jet axis near the jet edge 
and perpendicular towards the centre. Full line: data; dashed line: 
model. The profiles are blanked as in Fig. 1141 and the resolution 
is 0.75 arcsec. 



(v) The counter-jet polarization on-axis remains approx- 
imately constant where it can be measured (4 - 15 arcsec) 

(vi) Further than 4 arcsec from the nucleus, both the jet 
and counter-jet are polarized with a transverse apparent field 
on-axis. 

(vii) The apparent field is longitudinal along the edges of 
both jets wherever there is significant signal. 

(viii) There is a semicircular arc of enhanced polarized 
emission with a circumferential apparent field in the counter- 
jet, between 4 and 8 arcsec from the nucleus. 



4-3.5 Features that cannot be fitted well 
Our model cannot fit the following features: 

(i) The bright knot at the base of the main jet does not 
lie exactly on the jet axis and cannot be fitted precisely by 
our axisymmetric model. 

(ii) There is a second knot 1.9 arcsec from the nucleus, 
where the jet starts to expand rapidly (Fig. I12p and c). It 
would be possible to fit this feature by adjusting the emis- 
sivity profile to be much flatter locally, in which case the 
increase in path length through the jet would naturally gen- 
erate a maximum. In the absence of any constraint from the 
counter-jet, this complication cannot be justified. 

(iii) The arc of emission ~20 arcsec from the nucleus in 
the main jet, particularly noticeable in the polarized images, 
and similar to a prominent feature seen in 3C31 (LB), cannot 
be fitted. 
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Figure 16. Profiles of intrinsic parameters along the jets of 
B2 0326+39 as functions of z (not projected on the sky), (a) the 
shape of the jet edge (the scales are the same for both axes), 
(b) the velocity profile along the jet axis (solid line) and jet 
edge (dashed line), (c) the on-axis emissivity profile, converted 
to noB 1+a with no in m -3 and B in T. Solid line: model; dashed 
line: adiabatic approximation with the magnetic-field evolution 
expected from flux freezing. The normalization for the adiabatic 
profile is set so that it agrees with the free model at large distances 
from the nucleus. 



5 PHYSICAL PARAMETERS 

5.1 Model parameters and confidence limits 

The following section examines the structures of the jets 
implied by our models. All distances from the nucleus are 
given in linear units in the jet frame (i.e. not in projection 
on the sky). 
Table M 

gives the fitted parameters and error estimates, 
defined as in Section [3] for both sources. 
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Table 6. Fitted parameters and error estimates. 
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3.06 
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a The Symbol < means that any value smaller than the quoted maximum is allowed, 
k The Symbol > means that any value larger than the quoted minimum is allowed. 
Undetermined, as the counter-jet is not visible in this region. 



5.2 B2 0326+39 

5.2.1 Geometry and angle to the line of sight 

Our best-fitting model requires an angle to the line of sight 
of 64° ± 5°. The jet is initially well collimated, then flares 
between 3 and 8 kpc from the nucleus, reaching a maximum 
opening angle of ~20° (Fig. llGfa 'l. Further from the nucleus, 



it recollimates to a cone with an opening angle of only a few 
degrees. 



5.2.2 Velocity field 

The velocity profile along the jet axis is shown by the solid 
line in Fig. I16f b). The initial velocity of /3 ~ 0.8 is main- 
tained within 2 kpc of the nucleus. The jet then decelerates 
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uniformly out to ~7kpc from the nucleus, where the ve- 
locity reaches its asymptotic value. The data are consistent 
with a roughly constant velocity of /3 w 0.09 until the end 
of the modelled region. There are large fractional uncertain- 
ties, since variations in Doppler factor are slight, and we can 
only set an upper limit of /3 < 0.15 anywhere in this region. 
We discuss the velocity variations in the outer parts of the 
jet in more detail in Section fo. 21 

The transverse velocity structure is poorly constrained 
everywhere. In the high- velocity region close to the nucleus, 
the jet is narrow and poorly resolved transversely. At large 
z, the jet velocity is low, so differences in sidedness ratio 
between the centre and edge of the jet, which provide the 
principal constraint on transverse velocity variations, are un- 
detectably small for our preferred inclination. 



5. 2. 3 Emissivity 

The model on-axis profile of n B 1+a for B2 0326+39 (de- 
rived from the emissivity) is shown by the solid line in Fig 
lltif cl. In general, the profile flattens with increasing z: the 
fitted power-law indices are 6.1 (<1.0 kpc), 2.4 (1.0 - 2.7 kpc) 
and 1.3 (>3.0kpc). There are rapid changes in emissivity at 
1.0 kpc and 2.8 kpc. The former is explicitly modelled as a 
discontinuity, with the emissivity changing by a factor of 
1/g ss 40. The latter, fit by a power-law with a very steep 
index between 2.7 and 3 kpc, is also consistent with a dis- 
continuity. 

In the outer parts of the jet, the emissivity at the jet 
edge is «0.25 of its on-axis value with an error of a factor 
of «2. Elsewhere, the jet is poorly resolved transverse to its 
axis and the off-axis variation is poorly determined. 



5.2.4 Magnetic-field structure 

Fig. 1171 shows profiles of the radial, toroidal and longitudi- 
nal components of the magnetic field. As there is no varia- 
tion with streamline index, these curves apply to all stream- 
lines. The solid lines show our best-fitting model and the 
shaded areas the allowed errors. Close to the nucleus, the 
longitudinal component dominates. As the distance from the 
nucleus increases, this component drops, reaching a value 
consistent with zero by 8.5 kpc. The largest single compo- 
nent is toroidal from 3.5 kpc outwards, the fraction of ra- 
dial field remaining roughly constant everywhere. The field 
structure from 8.5 kpc outwards essentially consists of two- 
dimensional field sheets aligned perpendicular to the jet axis 
with a radial/toroidal field ratio «0.6. The transverse polar- 
ization profiles for the main and counter-jets (Fig. 1101 tell 
slightly different stories about the details of this field config- 
uration. In the main jet, the profile is essentially flat, and is 
marginally inconsistent with the model prediction of a higher 
degree of polarization on-axis. It would be slightly better fit 
by a model in which the toroidal and radial components are 
equal. The expected polarization for a non-relativistic jet 
with this field structure at an angle 8 to the line of sight 
is p = 0.43 for a = 0.55 llLainelll98fj) . This is slightly less 
than the observed value, and no transverse variation would 
be expected. In the counter-jet, the degree of polarization 
is higher on-axis (p ~ 0.6) and there are hints of a de- 
crease towards the edge. Both of these features are better 
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Figure 17. Profiles of: (a) radial, (b) toroidal and (c) longitu- 
dinal components of the magnetic field along the streamlines for 
B2 0326+39. Solid line: best fit model; shaded area: error on the 
best fit model derived from the limits in Table|S| dashed line: the 
profile expected on axis, in the absence of shear, if the magnetic 
field is frozen into the plasma. The latter curves are normalized 
to match the model at 6 kpc from the nucleus. 



fit by increasing the toroidal component. The optimization 
has compromised between the fits for the two jets to pro- 
duce the results shown in Fig. 1101 It is unclear whether the 
differences between the field component ratios in the main 
and counter-jets are real or a product of noisy data. 
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Figure 18. Profiles of intrinsic parameters along the jets of 
B2 1553+24 in the jet frame, (a) and (b) the shape of the jet 
edge. Panel (a) shows the jet with equal scales in both axes to 
give an idea of its real shape. Panel (b) has the width scale ex- 
panded by a factor of 6 to illustrate the variations in collimation. 
(c) the velocity profile along the jet axis (solid line) and jet edge 
(dashed line), (d) the on-axis profile of noB 1+a derived from the 
emissivity, with no in m — 3 and B in T. Solid line: model; dashed 
line, adiabatic approximation with the magnetic-field structure 
expected from flux freezing, normalized to match the free model 
at large distances. 

5.3 B2 1553+24 

5.3.1 Geometry and angle to the line of sight 

Our best-fitting model has the jet axis at 7.7° to the line 
of sight with an error of ±1.3°. This implies that the jet is 
has a very narrow 'pencil' shape, with a diameter of only 
^10 kpc at 200 kpc from the nucleus (Fig. 118b.). The jet is 
well collimated very close to the nucleus, flares to a maxi- 



mum opening angle of «6° between 15 and 35 kpc from the 
nucleus and then recollimates into a cone with an opening 
angle of «0.75° (FiglTBhV 

A near end-on orientation is qualitatively consistent 
with other measures (Section l2.3H and the source was indeed 
selected to be close to the line of sight. Nevertheless, the ex- 
tremely small value of 6 implies the existence of a large, 
missing population of side-on counterparts to B2 1553+24 
in the B2 sample. We address this topic in Section 1^31 



5.3.2 Velocity field 

Fig. HHf c) shows the model velocity profiles on-axis and at 
the jet edge (s = 1) for B2 1553+24. Within 20 kpc of the 
nucleus, the data are consistent with a constant value of 
(3 — 0.74, but the counter-jet is invisible in our images and 
we can only set a limit of f3 > 0.7. Between 20 and 30 kpc 
from the nucleus, the jet decelerates very rapidly and at 
large distances the velocity tends to an asymptotic value of 
P = 0.17 ±0.04. 

The variation in velocity transverse to the jet axis close 
to the nucleus is only constrained by the lack of a visible 
counter-jet, which indicates that the velocity at the jet edge 
must be (3 > 0.5 (> 0.7 times the on-axis value). Further 
from the nucleus, our best-fitting model has a ratio of edge 
to on-axis velocity «0.6, although the error is high. As in 
B2 0326+29, the velocities are too low to produce signifi- 
cant variations in Doppler factor across the jets unless the 
edge/on-axis ratio <0.2. 



5. 3. 3 Emissivity 

The on-axis profile of noB 1+a is shown by the solid line in 
Fig. I18l d1. Within ~5kpc of the nucleus the emissivity is 
constant, the dramatic brightening of the jet in this region 
being produced by the rapid expansion alone: no disconti- 
nuity is required. Further from the nucleus, the profile falls 
as a power-law, initially very steeply (index «8) but flatten- 
ing off as the distance increases. Between the start of the 
emissivity decrease and the end of the modelled region the 
emissivity drops by five orders of magnitude. In the regions 
of the jet where we have adequate transverse resolution, the 
emissivity falls to ~0.25 of its on-axis value at the jet edge 
with an error of roughly a factor of two. 



5.3.4 Magnetic-field structure 

Profiles of the rms field components along the jets are shown 
by the solid lines in Fig. ESI As in B2 0326+39, there is no 
variation with streamline index. The shaded areas show the 
variations allowed by the error estimates in Table |S| Within 
30 kpc of the nucleus, the field is mostly longitudinal; be- 
tween 30 and 40 kpc there is a rapid transition to a toroidally 
dominated structure which is maintained over the rest of the 
modelled region (although the error estimates do allow a sig- 
nificant proportion of longitudinal field to remain) . The best 
fit for the fractional radial component is close to a constant 
value of «0.2, but the data are also consistent with a com- 
plete absence of radial field. A pure toroidal field would fit 
the polarization data for z > 40 kpc. 
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Figure 19. Profiles of (a) radial, (b) toroidal and (c) longitu- 
dinal components of the magnetic field along the streamlines in 
B2 1553+24. Solid line: best-fitting model; shaded area: error on 
the model derived from the limits in Table 151 dashed line: profile 
expected on-axis, in the absence of shear if the magnetic field is 
frozen into the plasma. The latter curves are normalized to match 
the model values at 38 kpc from the nucleus. 



6 DISCUSSION 

6.1 Adiabatic models 

The simplest physical picture of the evolution of the emis- 
sivity along a jet assumes that the radiating particles only 
lose energy adiabatically (synchrotron and inverse-Compton 
losses being negligible by comparison) , that there are no dis- 
sipative processes such as particle acceleration or field-line 
reconnection, and that the magnetic field is convected pas- 
sively with the (laminar) flow. Analytical equations for the 



emissivity in this case wer e first derive d bv j Burchl il9T8f) : 
these were generalized by If3aum et alJ (Il997f) to describe 
a relativistic, decelerating flow with purely perpendicular 
or parallel field. We use the latter expressions, modified to 
calculate the m a gnetic -field evolution self-consistently as in 
lLaing fc Bridld (2004), to estimate the emissivity profiles 
expected for our models of jet shape and velocity. 

We refer to the fits derived in the previous sections and 
by LB as free models in order to distinguish them from adi- 
abatic models, following the terminology of lLaing fc Bridle! 
(2004). 



6.1.1 Magnetic- field profiles 

iBaum et al.l (119971) showed that the variations of the mag- 
netic field components in the quasi-one-dimensional approx- 
imation are: 



Br 

Bi 



in the absence of shear (x is again the jet radius). Given the 
field components at one point in the jet, we can then predict 
their evolution from the profiles of radius and velocity given 
in the previous section. 

The predicted evolution of the field components in 
B2 0326+39 and B2 1553+24 is shown by the dashed lines 
in Figs 1171 and 1191 respectively. The initial conditions have 
been chosen to match the free models at the distances from 
the nucleus given in the figure captions (note that the field 
is forced to be longitudinal at the nucleus). 

The qualitative evolution of the field components is con- 
sistent with the free models for both sources: a decrease in 
the longitudinal component is accompanied by an increase 
in the toroidal component. Rapid flaring and deceleration 
act in the same sense, leading to a large decrease in the rel- 
ative fraction of longitudinal field, as observed. There are 
significant quantitative discrepancies, however, particularly 
in the regions where the jets decelerate. In both sources, 
the transition from longitudinal to transverse field is pre- 
dicted to be less abrupt than is seen in the free models. 
This discrepancy is qualitatively consistent with the effects 
of velocity shear due to a transverse velocity gradient, which 
will act to increase the longitudinal component and to delay 
the onset of the longitudinal to transverse transition. This 
cannot be the whole story, however, as there are also large 
differences between the variations of the toroidal and ra- 
dial components. The former increases, as expected, but the 
latter stays roughly constant: in the quasi-one-dimensional 
approximation, their ratio sh ould remain constant . Similar 
pro blems occur in 3C31 (LB; lLaing fc Bridlejl2004tl . 

iKomissarov fc OvchinnikoT " F989lTl^9'ofi showed that 
the variation of the rms toroidal and longitudinal field com- 
ponents with distance could be very different from that 
predicted by the simple adiabatic appro ximation if the 
flow is turbulent , as is widely believed (e.g. lBicknelf"lll984t 
iDe Young"lll996l) . Given that the gross departures from the 
flux-freezing predictions occur where the jets are fast, or are 
decelerating rapidly, it seems likely that turbulence affects 
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both the magnitude and the configuration of the field in 
these regions. 



6.1.2 Emissivity profiles 

The emissivity function e can be written in terms of the 
magnetic field B, both as defined in Section f3,6l as: 

eoc(z 2 /?rr (1+2a/3) B 1+a 

jBaum et alJ Il997t lLaing fc Bridle! 12004) . B can be ex- 
pressed in terms of the parallel-field fraction / = (Bf) 1 ^ 2 / B 
and the radius x, velocity (3 and L orentz factor Y at som e 
starting location using equation 8 of lLaing fc Bridle! l|2u 



B oc 



1/2 



We can therefore predict the emissivity using our fitted 
jet width and velocity together with an estimate of the 
parallel-field fraction /. The resulting emissivity profiles for 
B2 0326+39 and B2 1553+24 are shown in Figs C^c) and 
I18f d~). respectively. The solid lines show the emissivity pro- 
files from our free model fits and the dashed lines the self- 
consistent adiabatic profiles, normalized to match at large 
z. 

In both objects the adiabatic models agree poorly with 
the free models where the jet velocities are high. The emis- 
sivity falls off much too rapidly in the adiabatic models, just 
as in 3C 31 (LB). This is not a surprise, for the following rea- 
sons: 

(i) Velocity shear (required by the free model) has been 
neglected. 

(ii) The magnetic-field evolution is more complicated 
than expected from simple flux-freezing in an axisymmet- 
ric laminar-flow model, even if shear is included (Sec- 
tion l6.1.T]l and turbulence may dominat e the field evolution 
jKomissarov fc Ovchinnikovlll989lll990l) . 

(iii) Where the jets are fast, we see complex, small-scale, 
non-axisymmetric structures, indicating that the flow is not 
laminar. 

(iv) In B2 1553+24, there is optical synchrotro n emission 
from the base of the main jet tParma et al.l l2003). implying 
continuing particle acceleration. 

Further from the nucleus, where the velocity has a low 
and roughly constant value, the adiabatic profile matches 
the free model very well in both objects. In B2 0326+39 this 
is within a region 4 to 12kpc (Fig. EE). In B2 1553+24, the 
region of low, constant velocity extends from 30 to 200 kpc 
and the adiabatic profile agrees very well with the free 
model over approximately 1.5 orders of magnitude in emis- 
sivity (Fig. I18H . This suggests that the outer jets in this 
source are modelled surprisingly well as constant-velocity, 
perpendicular-field, adiabatically-evolving flows. 



6.2 Sidedness profiles and reacceleration 

Our models, motivated by the gross features of the observed 
sidedness ratios, assume monotonic deceleration. There are 
theoretical reasons to expect jets to be reaccelerated by the 
pressure gradient of the external medium at large distances 



from the nucleus provided that the mass injection rate is not 
too large, for example if stel l ar mass loss dominates the mass 
inject ion llKomissarovlll994l : iBowman. Leahy fc Komissarovl 
1996). The most obvious effect of reacceleration is a small 
increase of sidedness ratio with distance from the nucleus. 
No such increase is obvious from the ridge-line profiles or 
images of sidedness ratio fFigslrjland ll3H . We have therefore 
averaged the sidedness ratios in rings of constant distance 
from the nucleus in order to improve the signal-to-noise ratio 

(Fig. HD . 

This reveals considerable sidedness structure in the low- 
velocity regions. Both sources show sidedness minima, at 
~14kpc from the nucleus in B2 0326+29 and at ~8kpc in 
B2 155 3+24. The minimum f o r B2 326+39 occurs roughly 
where IWorrall fc Birkinshawl fcOOOT) infe rred that the syn - 
chrotron minimum pressure in the jets feridle et al]ll99ll) 
becomes less than the pressure of the external medium. 

Fig. 1201 shows the average observed sidednesses as solid 
lines and our best-fitting models as dashed lines. The dot- 
ted lines show the profiles for models with the outer velocity 
parameters, j3o, modified to fit the minimum and maximum 
sections of the observed profiles in the outer region. The ve- 
locity ranges are (3 = 0.03 - 0.25 for B2 0326+39 and 0.13 
- 0.19 for B2 1553+24. [Note that the model sidedness pro- 
files increase slightly with distance from the nucleus because 
progressively larger areas of jet edge, which have lower sided- 
ness ratios than the centres, fall below the intensity blanking 
threshold] . 

If the changes in sidedness result from acceleration, then 
they should be associated with variations in the degree of po- 
larization. In B2 0326+39, the field structure can be roughly 
approximated by two-dimensional field sheets with equal ra- 
dial and toroidal components (Section 15. 2. 4H . The expected 
changes in polarization are then straightforward to calculate 
<Lainell98Cft . For an acceleration from (3 = 0.03 to /3 = 0.25 
between 16 and 18arcsec from the nucleus, as implied by a 
naive interpretation of Fig. l20f a'l. p should vary from 0.49 to 
0.64 in the main jet and from 0.45 to 0.38 in the counter-jet. 
These changes are not observed (FigsEt and|H). The velocity 
increase is far larger than is expected from the effects of any 
pressure gradient in the external medium and we also note 
that the fit of an accelerating adiabatic model to the emis- 
sivity profile would be significantly worse than that shown 
in Fig. EEI 

In B2 1553+24, it is more difficult to exclude accelera- 
tion as the cause of the increase in sidedness ratio between 
8 and 14arcsec from the nucleus (Fig. 120b ): the predicted 
variations in degree of polarization are below the fluctuation 
level in Fig. EM e ) an< I we cannot average across the jet to 
reduce the noise; the velocity increase required (/3o = 0.13 - 
0.19) is physically more reasonable than in B2 0326+39 and 
the adiabatic model fit to the emissivity is comparable in 
quality to that in Fig. 1181 

We conclude that there is no compelling evidence in 
favour of reacceleration of the jets in either source, and that 
it is most unlikely to be the sole cause of the increase in sid- 
edness ratio with distance from the nucleus in B2 0326+39. 
Nevertheless, it is expected theoretically, there are hints that 
it might occur, and it is consistent with the total intensity 
and polarization data in B2 1553+24. 
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Figure 20. Sidedness profiles of: (a) B2 0326+39 at 0.5 arcsec resolution and (b) B2 1553+24 at 0.75 arcsec resolution, averaging over 
all points at the same distance from the nucleus. Solid line: data; dashed line: best fit model; dotted lines: models with f3o modified to 
fit the minimum and maximum sidednesses in the lower-velocity region, as described in the text. 



6.3 The appearance of B2 1553+24 at other 
angles to the line of sight 

We infer that the angle to the line of sight for B2 1553+24 
is m8° . Given that the B2 sample from which it is drawn 
is selected at the low frequency of 408 MHz, the usual as- 
sumption is that the emission is isotropic and therefore that 
the distribution of jet orientations for its members should 
be random. The source should therefore have many (~100) 
counterparts of comparable total luminosity and larger 8, 
but the entire samp le of B2 radio galax ies with jets denned 
in Tables 1 and 2 of lParma et alJ <ll987l) has 43 members. A 
further concern is that the counterparts must be very large. 
The length of the main jet is at least 60 arcsec (Fig. [5^), cor- 
responding to 388 arcsec (340 kpc) at 8 — 60° , the median 
angle to the line of sight. This is compar able in size with th e 
longest jet in the sample, in NGC315 JWillis et alJJl98ll) . 
and far larger than the median («30 kpc: lParma et alll 987). 
There is cause for suspicion unless: (a) the side-on counter- 
parts of B2 1553+24 are not members of the B2 sample or 
(b) sources in that sample have linear sizes far larger than 
previously realised. 

In order to understand potential selection effects, we 
have computed the appearance of the model brightness dis- 
tribution for B2 1553+24 at various angles to the line of 
sight. We show the results for the median angle to the line 
of sight for an isotropic sample (8 — 60°) and for a source 
in the plane of the sky (8 = 90°) in Fig. HTJ The area cho- 
sen for the plots corresponds to the modelled region for the 
60° case. As well as being very long, the jets appear narrow, 
straight and (except for a small region around the nucleus) 
extremely faint. The total flux density from the inner 200 kpc 
of the model jets (measured along the jets and excluding the 



core) is constrained to be 44mJy for the 8 — 7.7° model but 
is 27mJy for 8 — 60°. There are three important selection 
effects: 

(i) The extended flux of the jets shows significant Doppler 
boosting at s mall 8. Given that there is little lobe emissi on in 
B2 1553+24 Jde Ruiter et alJll993tlCondon et al]ll99clh the 
jets probably dominate the flux at low radio frequencies, so 
end-on sources of this type will be selected preferentially in 
the B2 survey. The integral source count N(> S) oc S~ 3 ^ 2 
for the luminosities and redshifts in question, so the size 
of the parent sample is effectively increased by a factor of 
(44/27) 3/2 w 2 if isotropic lobe emission is ignored. 

(ii) Except for a small region around the nucleus, the 
model jets are faint. Further than 50 arcsec from the nu- 
cleus, the surface brightness is <0.16 mjy/beam area at 
8.4 GHz with a beam of 2.75 arcsec FWHM. This scales to 
0.47mJy/beam area at 1.4 GHz for a spectral index a = 0.6, 
comparable with the detection threshold f or typical VLA im- 
ages of sources in the sample presented bv lFanti et alJ il986L 
Il987l) . In any case, emission on scales > 120 arcsec would 
not have been imaged reliably in existing VLA observations 
of th e B2 sample (C configuration at 1.4 GHz; iFanti et alJ 
1987). The combination of these two effects makes it ex- 
tremely unlikely that the outer jets would have been de- 
tected. 

(iii) Sources of large angular size with flux densities close 
to the limit of the catalogue could have been missed by the 
original B2 survey, which measured peak rather than inte- 
grated flux density. We have estimated the magnitude of 
this effect by convolving the model brightness distribution 
for 8 = 60° with the beam of the B2 survey (3 x 10 arcmin 2 ). 
If the short axis of the beam is parallel to the jet axis, the 
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ratio of peak/total flux is 0.73; if the long axis of the beam is 
aligned, the ratio becomes 0.95. A few sources could there- 
fore have been missed (primarily those orientated E-W). 

We conclude that the side-on counterparts of B2 1553+24 
could have escaped identification, either because they were 
missed in the original survey or because their angular sizes 
have been greatly underestimated. Morganti & Parma (pri- 
vate communication) and Ledlow & Owen (2004; in prepara- 
tion) have made more sensitive observations of radio galax- 
ies in the B2 sample using the WSRT and the VLA in D 
configuration and have shown that a significant fraction of 
them have much longer radio jets than have previously been 
reported, extending many 100's of kpc or even further. We 
suggest that these form part of the missing population. 

6.4 Similarities and differences between 
B2 0326+39, B2 1553+24 and 3C 31 

In later papers, we will present a re-analysis of 3C 31 along 
with models of two further objects, all with the same func- 
tional descriptions of velocity, emissivity and magnetic field. 
This will allow a full comparison between them; here we 
make a few general observations about the similarities and 
differences between B2 0326+39, B2 1553+24 and 3C31. 
Notable similarities are: 

(i) All three objects can be modelled successfully by a 
decelerating jet model. 

(ii) The jets are all initially well collimated, then flare 
before recollimating and expanding conically. 

(iii) Close to the nucleus the jet velocity on-axis is con- 
sistent with a value of (3 ~ 0.8 in all cases. 

(iv) Deceleration occurs over a region of w 10 kpc to an 
outer velocity of 0.0 < f3 < 0.2. 

(v) The ratio of edge to on-axis velocity is consistent with 
a constant value of «0.7 everywhere, although in some places 
this is poorly constrained. 

(vi) The exponent of the emissivity index is similar in the 
outermost modelled regions of all the jets. 

(vii) At the jet edge the emissivity falls to ~0.25 of its 
on-axis value. 

(viii) The longitudinal field component dominates in the 
high-velocity regions close to the nucleus and the toroidal 
component in the outer parts, qualitatively but not quanti- 
tatively as expected from flux freezing. 

(ix) The variation of the radial field component is more 
complex than predicted by simple flux-freezing models. This 
effect cannot simply be due to velocity shear in a laminar, 
axisymmctric flow. 

(x) Close to the nucleus (in the high- velocity and deceler- 
ation regions), the variation of emissivity with distance from 
the nucleus is far less rapid than predicted by the adiabatic 
approximation. 

(xi) In the low-velocity outer regions, the emissivity vari- 
ation is reasonably well fitted by a quasi-one-dimensional 
adiabatic model. 

The three objects were expected to have a wide range of 
angles to the line of sight; this is confirmed by our modelling. 
In addition, there are some clear differences: 

(i) The conical outer region of 3C 31 has a large open- 
ing angle and is centred on the nucleus. The outer regions 



of the other two objects are closer to cylindrical, with very 
small opening angles and vertices far behind the nucleus. 
[Note that 3C 31 does recollimate outside the modelled re- 
gion (LB)]. 

(ii) The velocity in the outer region of 3C 31 was modelled 
by LB as decreasing monotonically with distance from the 
nucleus (as required by the variation of sidedness ratio). In 
the B2 sources, it has an approximately constant asymptotic 
value. 

(iii) The close connection between the jet geometry and 
the forms of the velocity and emissivity profiles inferred for 
3C 31 is not general. 

(iv) The velocity in B2 1553+24 remains at f3 ^ 0.8 un- 
til «20kpc from the nucleus (cf. ~2kpc in the other two 
objects). 

(v) There is no need for a discontinuity in emissivity in 
B2 1553+24: the brightening of the jet is consistent with 
expansion at constant emissivity. 

(vi) Although the biggest single field component at 
large distances is always toroidal, there are significant dif- 
ferences in the details of the field-component evolution: 
3C 31 has a mixture of longitudinal and toroidal compo- 
nents; B2 0326+39 has toroidal and radial components and 
B2 1553+24 has almost pure toroidal field. 

(vii) In 3C31, the radial field component in the flaring 
region increases towards the edge of the jet. There is no 
evidence for this effect in the B2 sources, although it cannot 
be excluded in B2 0326+39, where the signal-to-noise ratio 
is low. 

(viii) The adiabatic approximation describes the emissiv- 
ity evolution in B2, 1553+24 over a much larger fraction of 
the jets than in the other two sources. This may be because 
it applies most accurately in the low- velocity outer region, 
which is relatively longer in this source. 

7 SUMMARY AND FURTHER WORK 
7.1 Summary 

7.1.1 General 

We have shown that the total and linearly-polarized syn- 
chrotron emission from the jets of the FRI radio galaxies 
B2 0326+39 and B2 1553+24 can be modelled successfully 
on the assumption that they are axisymmetric, intrinsically 
symmetrical, relativistic, decelerating flows. The models are 
based on those developed by LB for 3C 31, with modifica- 
tions to the functional forms assumed for geometry, velocity, 
emissivity and magnetic field. In particular, we found that 
the variations of intrinsic parameters were not directly cou- 
pled to the geometry of the jets as defined by their outer 
isophotes. 

7.1.2 Geometry 

In the flaring region close to the nucleus, the jet radius x 
can be fitted by a polynomial of the form x = a^z 2 + 03 z 3 , 
where z is the distance from the nucleus. In the outer region 
at larger distances, the jets are conical, with opening angles 
sufficiently small that they are almost cylindrical. The in- 
ferred angles to the line of sight are 64° ± 5° for B2 0326+39 
and 7.7° ± 1.3° for B2 1553+24. 
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Figure 21. The brightness distribution of the model for the jets in B2 1553+24 at angles to the line of sight of: (a) 9 = 60° and (b) 
9 = 90°. The plots cover 394 arcsec (346 kpc) in projection on the sky. This corresponds to a distance from the nucleus along the jet axis 
of 200 kpc in both directions for 9 = 60°, i.e. to the size of the modelled region. The contour levels are 1, 2, 4, 8, 16, 32, 64, 128, 256, 
512 X 20/jJy/bcam and the model images have been convolved with a 2.75-arcsec Gaussian beam. 



7.1.3 Velocity 

Close to the nucleus the jets have high and approximately 
constant velocities /3 ~ 0.8. They decelerate rapidly over dis- 
tances of < 10 kpc to a constant asymptotic outer velocities 
(3 < 0.2. The velocities at the jet edges are consistent with 
«0.7 of their on-axis values everywhere. 

7.1.4 Emissivity 

The emissivity profiles along the jets can be divided into dis- 
crete regions, each with a power-law variation. The power- 
laws tend to be natter further from the nucleus, although 
regions close to the nucleus (where knots of emission are 
seen) are modelled as bright plateaux with rapid decreases 
in emissivity at their ends. In the low-velocity regions after 
the jets have decelerated, the emissivity profiles are consis- 
tent with a simple, quasi-one-dimensional adiabatic model; 
closer to the nucleus, the adiabatic approximation predicts 
too steep a decline of emissivity with distance. The fractional 
emissivity at the edge of the jets is ~0.25 everywhere. 

7.1.5 Magnetic Field 

Although the observed polarization structures of the two 
objects appear radically different, the main features of the 
magnetic-field evolution are similar in both cases: the prin- 
cipal differences result from projection. In both objects, the 
field is dominated by the longitudinal component in the 
region close to the nucleus; further out, the largest single 
component is toroidal. The radial field components are ~0.4 
and Ri0.2 of the total in B2 0326+39 and B2 1553+24, re- 
spectively. No evidence for transverse structure in the field- 
component ratios is seen in either source. Evolution from 
longitudinal to transverse field is qualitatively consistent 
with flux-freezing in an expanding, decelerating jet, but 
there are complications. First, the ratio of radial to toroidal 
field (which should remain roughly constant even if there 
is velocity shear) varies significantly along the jets. Second, 
the details of the component evolution are quantitatively 
inconsistent with the flux-freezing picture. 

7.2 Further work 

We have modelled 3C 31 with the modified functional de- 
scriptions used here in order to allow a quantitative com- 



parison with other objects. As expected from the variation 
of sidedness ratio with distance from the nucleus (LB), the 
jet velocity in 3C 31 does not reach an asymptotic value 
in the modelled region. We have shown, however, that use 
of the new velocity law allows as good a fit to the ob- 
servations as that achieved by LB and that their conclu- 
sions are not significantly modified. We have also success- 
fully appl ied the model to t he well-know n radio galaxies 
NGC 315 JVenturi et alJll993tl and 3C 296 faarcastle et alJ 
Il997fl . These results will be presented in future papers. 

When suitable X-ray observations are available, we 
will apply the conser vation-law approach developed by 
lLaing fc Bridlel \ 2002lJl to infer the energy and momentum 
fluxes of the modelled jets and their variations of pressure, 
density and entrainment rate with distance from the nucleus. 
Our eventual aim is to replace the empirical descriptions of 
internal quantities with more physically derived ones. Ini- 
tially, we will apply the more sophis ticated adiabatic model 
developed bv lLaing fc Bridlel i2004t) to the outer regions of 
the jets and in particular to B2 1553+24, where the simple 
analysis of Section [6.11 suggests that the adiabatic approxi- 
mation holds over a large range of distances. We then plan 
to incorporate the effects of synchrotron and inverse Comp- 
ton losses in a self-consistent way, to model the observed 
emission from radio to X-ray wavelengths and to quantify 
the particle-acceleration processes in FRI radio jets. 
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